Lithological tomography with the correlated pseudo-marginal method
نویسندگان
چکیده
We consider lithological tomography in which the posterior distribution of (hydro)geological parameters interest is inferred from geophysical data by treating intermediate properties as latent variables. In such a variable model, one needs to estimate intractable likelihood given data. The pseudo-marginal method an adaptation Metropolis-Hastings algorithm unbiased approximation this obtained Monte Carlo averaging over samples from, setting, noisy petrophysical relationship linking and properties. To make practical data-rich settings with low noise levels, we demonstrate that sampling must rely on importance distributions well approximate scatter around sampled parameter field. achieve suitable acceptance rate, both (1) correlated method, correlates used proposed current states Markov chain, (2) model proposal scheme preserves prior distribution. As synthetic test example, infer porosity fields using crosshole ground-penetrating radar (GPR) first-arrival travel times. use (50x50)-dimensional pixel-based parameterization multi-Gaussian field known statistical parameters, resulting space high dimension. our prior-preserving outperforms state-of-the-art methods linear non-linear greatly enhancing exploration.
منابع مشابه
Pseudo-Marginal Slice Sampling
Markov chain Monte Carlo (MCMC) methods asymptotically sample from complex probability distributions. The pseudo-marginal MCMC framework only requires an unbiased estimator of the unnormalized probability distribution function to construct a Markov chain. However, the resulting chains are harder to tune to a target distribution than conventional MCMC, and the types of updates available are limi...
متن کاملLearning Gaussian Graphical Models With Fractional Marginal Pseudo-likelihood
We propose a Bayesian approximate inference method for learning the dependence structure of a Gaussian graphical model. Using pseudo-likelihood, we derive an analytical expression to approximate the marginal likelihood for an arbitrary graph structure without invoking any assumptions about decomposability. The majority of the existing methods for learning Gaussian graphical models are either re...
متن کاملPseudo-acoustic wavefield tomography with model constraints
Anisotropic wavefield tomography seeks to recover anisotropic parameters based on seismic wavefields. In terms of computational efficiency and degrees of freedom of the model parameters, the acoustic approximation to the anisotropic elastodynamics is more suitable than the full anisotropic elastic equations to the anisotropic tomography problem. We implement the wavefield tomography in VTI medi...
متن کاملMarginal methods for correlated binary data with misclassified responses
Misclassification is a longstanding concern in medical research. Although there has been much research concerning error-prone covariates, relatively little work has been directed to problems with response variables subject to error. In this paper we focus on misclassification in clustered or longitudinal outcomes. We propose marginal analysis methods to handle binary responses which are subject...
متن کاملذخیره در منابع من
با ذخیره ی این منبع در منابع من، دسترسی به آن را برای استفاده های بعدی آسان تر کنید
ژورنال
عنوان ژورنال: Geophysical Journal International
سال: 2021
ISSN: ['1365-246X', '0956-540X']
DOI: https://doi.org/10.1093/gji/ggab381